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XMM- Newton observation of the bursting 
pulsar GRO J1744-28 in quiescence 
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Abstract. The XMM-Newton X-ray observatory performed a pointed observation of the bursting pulsar GRO 
J1744-28 in April 2001 for about 10 ks during a program devoted to the scan of the Galactic center region. 
After the discovery of this source by BATSE in December 1995 during a very active bursting phase, it has been in 
quiescence since April 1997. We present here the first detection of GRO J 1744-28 in its quiescent state at a position 
which is consistent with previous high-energy positions (but not consistent with a proposed IR counterpart). The 
observed luminosity of the source in quiescence is about 6 orders of magnitude weaker than the luminosity in 
outburst. 
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1. Introduction 

The bursting pulsar GRO J 1744-28 is a unique source 
that presents both the characteristics of a type II X- 
ray burster and a X-ray pulsar. It was discovered near 



the GalaeLie center by Lhe BuisL and Transient Sumee 
Experiment (DATSE) uu buaid Lhe GumpLuii Gaiiima- 

Ray Observatory (CGRO) on 1995 December 2 ( Fishmar] 
|ct al. 1995 ; Kouveliotou ct al. 1996bj ), when it was ex- 
periencing a very active bursting phase with a rate of 
about 20 hard X-ray bursts per hour. After one day, 
the source entered a regime of hourly bursting which 
lasted until 1996 April 26. Many of these bursts were 
also detected by the Proportional Counter Array (PCA) 
instrument on board the Rossi X-ray Timing Explorer 



GRO J 1744-28 is a 467 ms X-ray pulsar in a binary system 
of orbital period 11.8 days. The pulsar period determina- 
tion was confirmed with ASCA ( [Dotani et al. 1996a| ]b|). 
It has also been found that the X-ray bursts were show- 
ing pulsations with the same 467 ms period (Kouveliotou 



ct al. 1996c) and had similar hard X-ray / soft gamma-ray 



spectra than the persistent emission (Briggs et al. 1996; 
[jtrickman et al. 1996] ), so that the association is firmly 
established. 

The initial source localization was a 6° radius error circle 



(Fishman et al. 1995), which was rapidly refined to a 24' 



(RXTEj-( Swank 1990 r Giles et al. t ggt ) and Lhc Oriented 



wide annulus by the Interplanetary Network using time- 
delay measureme nts between CGRO and other satellites, 
mainly Ulysses ( Burley et al. 1995 ). Obser vations with 
ASCA provided a ~ Y error circle. ROSAT ( Kouveliotov, 



ScmtilMioii Specb ronieLer Experiment. (O SSE) InsU'iimeiil 
on board CGRO (3trickman et al. 1996). Then the burst 
rate decreased rapidly and the bursting activity ceased for 
the next 7 months. On 1996 December 1, a second burst- 
ing phase started with a similar burst rate to that seen in 
the first outburst flWoods et al. 1999j ). It lasted until 1997 
April 7. The source was no longer detected by BATSE or 
RXTE after this date. More than 5800 bursts were ob- 
served by BATSE in total. 

During the first outburst, persistent hard X-ray emission 
f rom GRO J1744-28 was detected by BA TSE and OSSE 
( Paciesas et al. 1996 ). Finger et al. (1996 ) measured that 



et al. 1996a ; Augusteijn et al. 1997 ) improved the source 
localization to R.A. 17 h 44 m 33M and decl. -28°44'29" 
(J2000) with a ~ 10" error circle. An optical / near- 
infrared variable object was observed in 1996 at the pe- 



riphery of this ROSAT error circle (Augusteijn et al. 1997; 
pole ct al. 1997 ) and proposed as a possible counterpart. 
More recently, Hurley et al. (2000) derived a more pre- 
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cise IPN localization of GRO J1744-28 with a 3cr error 
ellipse of 532 arcsec 2 fully contained in the ASCA error 
circle and partially overlapping with the ROSAT error cir- 
cle. The optical / near-infrared variable source lies outside 
the new IPN 3<r ellipse, which makes it unlikely to be as- 
sociated with GRO J1744-28. 

Finger ct al. (199(f ) determined the mass function of the 
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Inst. 


RA 


DEC 


stat. 


error 


error 








error 


(RA) 


(DEC) 


PN 


17 h 43 m 51 s .3 


-28° 46' 40" 


~ 0.6" 


-0.3" 


+1.8" 


M0S1 


17 h 43 m 51 B .4 


-28°46'38" 


~ 0.7" 


-1.2" 


+0.1" 


MOS2 


17 h 41 m 51 s .4 


-28°46'37" 


~ 0.7" 


-1.7" 


-0.8" 



Table 1. Position of the reference star TYC 6840-38-1 
as detected by each of the EPIC instruments. The first 
three columns (right ascension, declination and statisti- 
cal error at la level) are the result of the standard SAS 
source detection procedure. The last two columns show 
the measured offsets with respect to the Tycho-2 cata- 
logue position R.A. 17 h 43 m 51 s .3 and decl. -28°46'38". 



Inst. 


RA 


DEC 


stat. 
error 


counts 


PN 


17 h 44 m 33 s .l 


-28°44'26" 


~ 2.4" 


22.4 + 6.1 


MOS1 


17 h 44 m 33 s .2 


-28°44'27" 


~ 1.9" 


27.4 + 6.0 


MOS2 


17 h 44 m 33 s .l 


-28°44'16" 


~ 6.0" 


18.8 + 7.4 


EPIC 


17 h 44 m 33 s .2 


-28°44'26" 


~ 1.4" 


59.2 ± 10.0 



Table 2. Position of GRO J 1744-28 as detected by each of 
the EPIC instruments. The right ascension (RA), declina- 
tion (DEC), statistical error at la level and counts number 
are the result of the SAS standard source detection proce- 
dure. The last line (EPIC) corresponds to the image made 
of the sum of the PN, MOS1 and MOS2 images. 



binary system to be 1.36 10~ 4 M©. This supports that 
GRO J 1744-28 is very likely a Low-Mass X-ray Binary 
(LMXB), where the central neutron star accretes mat- 
ter through Roche lobe overflow from an evolved low- 



mass stellar companion (Daumcric ct al. 1996; Sturner & 



Dcrmc " 1996). This is however the only known example of 
such a system showing both the characteristics of a X-ray 
pulsar (the persistent pulsed emission) and of a type II 
X-ray burster. The relative fluences in the bursts and in 
the persistent emission probably rule out a thermonuclear 
origin of the bursts, which are more likely of type II. Such 
bursts are accretion powered and are possibly triggered 
by an instability in the accretion flow, the nature of which 



is not clearly identified (Kouveliotou et al. 1996c; Lewin 



et al. 1996; Cannizzo 1996). This is supported by several 
similarities between the bursts of GRO J1744-28 and those 
of the Rapid Burster MXB 1730-355 QLcwin ct al. 1996| ), 

1. e the first and unique object where type II bursts were 
discovered. 

In this paper, we report new observations of GRO J1744- 
28 in the X-ray band, made with XMM-Newton in April 
2001 as it was in its quiescent state. 

2. Observations and data analysis 

GRO J1744-28 was observed with the European Photon 
I maging Camera (E PIC), both with the MOS cameras 
( [Turner et al. 20014) o n 2001 April 4 at UT 13:14:27 and 
with the PN camera ( Struder et al. 2001 ) on the same 
day starting at UT 14:27:23. The MOS and PN cameras 
were operating in full window mode with the medium fil- 
ter. Data reduction was performed using the XMM SAS 
(Science Analysis Software) and the standard XSPEC 
packages. 

The time histogram of all events above 10 keV showed 
that the background has only small variations, without 
any very intense flare, even if the averaged flux is close to 
the alert level as defined in the standard SAS flare rejec- 
tion procedure. For this reason we applied no time filtering 
so that the total duration has been kept to 9199 s for the 
MOS 1 & 2 cameras and 4699 s for the PN camera. A pre- 
liminary visual inspection of the MOS 1,2 and PN images 
showed that a faint source was present at the position of 
GRO J1744-28 (with no significant flux below 1 keV and 



above 6 keV). The absence of flux below 1 keV can be 
explained by the strong absorption towards the Galactic 
center. Therefore all the further analysis has been made in 
the 1-6 keV band. Also, it could be noticed that for some 
unidentified reasons, the background in the MOS2 image 
appears to be higher than in the MOS1 image. 
In this energy band, we generated an image for each of 
the three EPIC instruments, keeping only the valid X-ray 
events using the canonical range of pattern specified by 
the SAS (0-12 for the MOS cameras and 0-4 for the PN 
camera). We applied to these three images the standard 
source detection procedure with the following SAS tasks : 
(i) eexpmap and emask to generate the exposure map and 
the detection mask for each image; (ii) eboxdetect (in local 
mode) to generate an input list of source positions; (iii) es- 
plinemap to generate the corresponding background map; 
(iv) eboxdetect (in map mode) to generate a new list of 
sources with a better detection sensitivity using the back- 
ground map; (v) emldetect to perform a maximum like- 
lihood PSF fit to the count distribution of each source. 
This procedure was also applied to an image made of the 
sum of all X-ray events of the three instruments in the 1-6 
keV band. 

The standard SAS source detection procedure that we ap- 
plied to the data found one weak source within the ROSAT 
~ 10" error circle of GRO J 1744-28. This source was also 
present in the standard SAS pipeline output. As can be 
seen in Table |[ it is a detection at the ~ 4a level for the 
PN and MOS1 cameras and a detection at the 2.5a level 
only for the MOS2 camera. The number of sources at the 
3a level in the total field of view of the MOS cameras 
(radius 15') is about 100 so that the probability to find 
such a source by chance in an area comparable with the 
ROSAT error box (radius - 10") is about 1 % . Then we 
identify this weak source with GRO J1744-28. 
The position error obtained with the standard SAS source 
detection procedure can be dominated by the errors of 
the absolute astrometry (maximum ~ 4"). To estimate 
these systematic errors, we proceeded to identify strong 
sources in our field of view. Unfortunately only one is 
strong enough to be used in this way. This bright source 
which was detected by each instrument has been identified 



Daigne et al.: XMM observation of GRO J 1744-28 in quiescence 



3 



based on positional coincidence - as the star TYC 6840- 



Spectrum of GRO J1744-28 



38-1 in the Tycho-2 Catalogue (Hog et al. 2000) which con- 
tains the position, the proper motion and the two-colour 
photometric data for the 2.5 million brightest stars in the 
sky. Its catalogue position is R.A. 17 h 43 m 51. s 3 and decl. 
— 28°46'38" (with a negligible proper motion for our pur- 
pose). Table [l] shows the position of this star obtained 
with each instrument. They are all consistent with the 
catalogue position, taking into account the statistical er- 
ror given by the emldetect procedure (~ 0.7") and the 
quoted accuracy of the XMM-Newton coordinates: the 
MOS1 and MOS2 positions are distant of ~ 1" from each 
other whereas the PN position is shifted of ~ 3" with re- 
spect to the MOS2 position. This is in agreement with 
the negligible offset found in the MOS 1&2 roll angle and 
t he mean offset of ~ — 0.30° found in the PN roll angle 
(fTedds & Watson 200l|). As the star is located at ~ 10' 



of the center of the field of view, the expected shift in the 
PN position is indeed ~ 3". However the same shift at the 
position of GRO J1744-28 is negligible because the source 
is at the center of the field of view. The nice agreement 
bewteen the position of the star obtained with the three 
EPIC instruments and the catalogue position shows that 
for our observation, the systematic error in the position is 
~ 1". 

The position of GRO J1744-28 is the same in the three 
instruments taking into account the statistical error (see 
Table |). The best accuracy is obtained when summing 
the events of the three instruments. The derived position 
is R.A. 17' l 44 m 33 s .2 and decl. -28°44'25" with an error 
circle of about 4" at the 3tr level (summing in quadrature a 
statistical error of ~ 4" and a systematic error ~ 1"). The 
closest other source is located at ~ 25" of GRO J1744-28, 
well outside the ROSAT error circle, and therefore can be 
easily rejected. 

To construct the spectral distribution of the photons, we 
have decided to sum the two MOS cameras events to in- 
crease the statistics, which is possible because the response 
matrix are similar. We extracted from the PN and MOS 
1+2 images the X-ray events within a radius of 15" of the 
GRO J1744-28 position. We grouped the photons in 5 en- 
ergy bands : 1 - 2 keV, 2-3 keV, 3-4 keV, 4-5 keV and 
5 — 6 keV. Background was estimated using offset regions 
in the same central CCD. After substracting it, we ob- 
tained a detectable flux only in the 1 — 5 keV band. The re- 
sulting spectrum was obtained with XSPEC and is plotted 
in Fig. [j]. The source count rate in the 1-6 keV energy band 
is (4.75 ±1.32) 10 -3 counts/s for the PN camera and 
(3.40 ± 0.76) 10~ 4 counts/s for the MOS 1+2. Fixing the 
hydrogen column density to the value of 5.1 x 10 22 cm~ 2 
(Dotani et al. 1996a), we tried to fit the spectrum of GRO 
J1744-28 either with a power-law or with a blackbody dis- 
tribution. We found that acceptable fits of the observed 
spectrum with a power-law could be obtained for a pho- 
ton index in the range ~ 2 — 5. Acceptable fits with a 
blackbody distribution could be obtained as well for a 
temperature kT in the range ~ 0.4 — 1. keV. The cor- 
responding unabsorbed energy flux in the 1 — 5 keV is 



energy (keV) 



energy (keV) 

Fig. 1. Spectrum of GRO J1744-28 : the flux is plotted as 
a function of energy for the PN camera (top) and the sum 
of the MOS1 and MOS2 cameras (bottom). The solid line 
corresponds to the fit of the data by an absorbed black- 
body spectrum with a column density 5.1 10 -22 cm~ 2 and 
a temperature kT — 0.6 keV. 



about 2 - 
about 1 



5 10" 
- 2 10 



1 S — 2 

erg. cm .s 



-13 



erg. cm 



1 in the power-law case and 
_1 in the blackbody case. 



3. Results and discussion 

3.1. Source position 

Fig. H shows the XMM-Newton error circle we have de- 
rived in the last section. It is entirely included in the in- 



tersection of the ROSAT error circle (Augusteijn et al 



1997) and the IPN error ellipse (Hurley et al. 200C). This 



strengthens the identification of the source observed by 
XMM-Newton with GRO J1744-28, despite the fact that 
this source is too weak to allow the search for pulsations 
(such pulsations were also impossible to study with the 
ROSAT observation). The association of the proposed op- 
tical / near-infrared counterpart (Augusteijn et al. 1997; 
pole et al. 1997| ) with GRO J1744-28 is rejected without 
any ambiguity with the more precise XMM-Newton er- 
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ROSAT 



17:44:34 17:44:33 17:44:32 17:44:31 
RIGHT ASCENSION 



17:44:30 



Fig. 2. An ~ 1' x 1' square region containing the approx- 
imate 3tr IPN error ellipse, the ~ 10" radius ROSAT er- 
ror circle and the 3<r XMM error circle (~ 4" radius). 
The two positions of the proposed optical / near-infrared 
counterpart (Augusteijn et al. 1997; Cole et al. 1997) are 
respectively marked with a star and a cross. 



ror circle we have obtained in the previous section. It is 
then important to carry out deeper searches for the real IR 
counterpart of GRO J1744-28 using the better determined 
XMM-Newton position. 

3.2. Source luminosity in quiescence 

Despite a very low statistics, we tried in the last section 
to estimate the spectrum of GRO J1744-28 in quiescence. 
This spectrum appears to decrease very fast at high en- 
ergy, as it is indicated by the photon index found when 
fitting with a power-law. The unabsorbed luminosity in 
the 1-5 keV band is 



(1 



6) 10 33 erg/s 



(1) 



where d\o is the distance of the source in unit of 10 kpc. 
The quiescent luminosity is then about five orders of mag- 
nitude smaller than the outburst luminosity 



L 



10 38 d^ 



erg/s 



(2) 



and is compatible with the typical quiescent luminosity of 
a neutron star transient ( Asai et al. 1998| ). On the other 
hand, if we assume a blackbody distribution, the range 
in temperature we obtain (0.4-1. keV) is slightly outside 
the typical range of neutron star tran sients in quiescence 
(kT ~ 0.2 - 0.3 keV; |Asai et al. (1998[) ). The temperature 
we find corresponds to an emitting area 



.4 



L 
^f 1 



(10 9 - 10 11 ) d^ cm 2 



(3) 



Such an area corresponds to a sphere of radius 10 4 — 
— 10 5 cm which is well below the size of a neutron star. 
Then the observed luminosity cannot be associated with 
the thermal emission of the total surface of a non accreting 
neutron star. Many other origins can be proposed, such as 
the emission due to the accretion on the magnetic poles 
or the emission from the accretion onto the neutron star 



magnetosphere (Campana et al. 1998). However the large 



uncertainty we get on the spectral distribution of the X- 
ray photons makes very difficult any detailed comparison 
with a particular model. 
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